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Accelerating polarized proton beams and storing them for many turns can lead to a loss of 
polarization when accelerating through energies where a spin rotation frequency is in resonance 
with orbit oscillation frequencies. First-order resonance effects can be avoided by installing Siberian 
Snakes in the ring, devices which rotate the spin by 180° around the snake axis while not changing 
the beam's orbit significantly. For large rings, several Siberian Snakes are required. 

Here a criterion will be derived that allows to find an optimal choice of the snake axes. Rings 
with super-period four are analyzed in detail, and the HERA proton ring is used as an example 
for approximate four-fold symmetry. The proposed arrangement of Siberian Snakes matches their 
effects so that all spin-orbit coupling integrals vanish at all energies and therefore there is no first- 
order spin-orbit coupling at all for this choice, which 1 call snakes matching. It will be shown that 
in general at least eight Siberian Snakes are needed and that there are exactly four possibilities 
to arrange their axes. When the betatron phase advance between snakes is chosen suitably, four 
Siberian Snakes can be sufficient. 

Since the spin motion depends on a particle's trajectory, protons at different phase space positions 
generally have different spin directions. The time averaged polarization at each phase space point is 
parallel to the invariant spin field, and the spread of this spin field limits the polarization that can 
be stored. By the here presented choice of Siberian Snakes this limit is completely eliminated up 
to first order in the transverse phase space coordinates. The invariant spin field and the amplitude 
dependent spin tune are also computed without linearization to show the advantages of this choice 
of snakes. Ultimately, the goal of snake matching is to reduce the loss of polarization during the 
acceleration of the beam. To show that favorable choice of snakes have been found, polarized protons 
are tracked for part of HERA-p's acceleration cycle which shows that polarization is preserved best 
for the here proposed arrangement of Siberian Snakes. 



I. INTRODUCTION 

The design orbit spin direction uq at some azimuth 
9 = 2irl I L along a storage ring of length L describes the 
spin direction for the design orbit which is periodic from 
turn to turn, where / is the pathlength along the design 
orbit. The design orbit spin tune v$ describes the number 
of times a spin has rotated around uq{6) during one turn 
around the ring. 

Spin-orbit resonances in high energy accelerators arise 
when the electro-magnetic fields on synchro-betatron tra- 
jectories cause disturbances of the spin's motion which 
build up coherently from turn to turn. 

In a flat ring, for instance, an initially vertical spin of 
a particle traveling on the design orbit remains vertical 
during particle motion. On a vertical betatron trajectory 
the particle traverses horizontal fields in quadrupoles and 
the spin no longer remains vertical. This disturbance of 
spin motion due to the betatron motion is described by 
the spin-orbit coupling integrals Q] 

/± = f n kyS-^yUe , (1) 

Jo 

where k y — with the quadrupole strength k and 

the vertical beta function (3 V , and '5 is the phase advance 
of the spin rotation around the vertical and $ a is the 
vertical betatron phase. For ease of notation a constant 
factor which is sometimes multiplied to this integral is 
not indicated here. 



The I tend to be especially big when the spin distur- 
bance in every FODO cell of a ring adds up coherently 0] • 

In an approximation that is of first order in transverse 
phase space coordinates, these integrals yield the follow- 
ing important information: If all spin-orbit coupling inte- 
grals vanish, all initially vertical spins are again vertical 
after one turn, although they have traveled along differ- 
ent betatron trajectories. The ring is then called spin 
matched or spin transparent. 

In this paper it will be analyzed how Siberian Snakes 
can be used to make a storage ring spin transparent at 
all energies. 

When a beam is polarized, it has a polarization di- 
rection f(z,j) at each phase space position z at some 
azimuth 9 of the ring, where j indicates how often the 
bunch has traveled around the ring. Since each initial 
phase space position leads to a different particle trajec- 
tory with different spin motion, spin fields in general 
change from turn to turn. A special spin field that does 
not change from turn to turn is called the invariant spin 
field n{z) @. 

When a beam is polarized according to this spin field, 
f(z, j) = n(z), then the beam polarization is given by 
the phase space average Pu m =< n{z) >g. In any case, 
even when a particle at z{j) after its jth turns has a clas- 
sical spin vector S with a direction which is not parallel 
to n(z(j)), its projection onto the invariant spin field is 
constant, i.e. f(z(j),j) ■ n(z(j)) does not change with j. 
This is possible since its spin S appears to rotate around 
n(z(j)). The time averaged polarization at the phase 
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space point z is therefore parallel to n(z) for any spin 
field. The maximum time averaged polarization at that 
point is thus given when the beam is initially polarized 
parallel to n(z). The average Pum is therefore called the 
limiting polarization 0. 0- Q . 

To illustrate one of the benefits of the proposed choice 
of Siberian Snakes it will be shown that they increase 

Plira • 

When an invariant spin field is found, the rotation of 
S around it can be described by an amplitude dependent 
spin tune v(J) that depends on the orbital amplitudes J 
but not on the orbital phase variables. This is important 
since v(J) does not change from turn to turn and can 
therefore be used to describe long term coherence with 
the frequencies of orbital motion. When this amplitude 
dependent spin tune is in resonance with the orbital tunes 
Qj,i.e. v{J) = k + J^jkjQj f° r integers kj, the spin 
motion can be strongly disturbed. 

As another illustration of the benefits of the proposed 
choice of Siberian Snakes it will be shown that they in- 
crease the orbital amplitudes for which low order reso- 
nance conditions occur. 

It can be shown that the projection of the polarization 
onto the invariant spin field J$ = S ■ n(z) is not only 
invariant in a storage ring, but that it is an adiabatic 
invariant when parameters of the accelerator, like the 
beam's energy, change. This means that a beam that 
is initially polarized parallel to the invariant spin field 
n(z,Ei) and thus has the average polarization Pu m (Ei) 
at its initial energy Ei will still be polarized parallel to 
the invariant spin field n(z, Ef) after the beam has been 
accelerated to a final energy Ef, and its average polar- 
ization is Pn m (Ef). This final polarization can be very 
large, even though it may have been very small at some 
intermediate energies. All this is only true if the change 
of energy is performed adiabatically slowly. Usually, the 
acceleration cannot be performed adiabatically slowly at 
all energies and polarization is lost, i.e. the final polar- 
ization is smaller than Pu m (Ef). 

The most convincing illustration of the benefits of the 
proposed choice of Siberian Snakes is finally that they 
increase the polarization that is retained after accelera- 
tion. 

Although it has been straightforward to define n(z), 
it is not easy to calculate this spin field in general and 
much effort has been spent on this topic, mostly for 
electrons at energies up to 46 GeV. All algorithms de- 
veloped before the polarized proton project at HERA-p 
rely on perturbation methods at some stage, and either 
do not go to high enough order [H, |t| or have problems 
with convergence at high order and high proton energies 
[H [H EJ- The algorithms developed for the HERA-p 
project P.l5l lT^.IT3 . ll5| made the here presented analysis 
possible. 



II. OPTIMAL CHOICES OF SIBERIAN SNAKES 

Siberian Snakes are indispensable if polarized proton 
beams are to be accelerated in a high-energy synchrotron 
such as HERA-p. This has the following reasons: 

1. Siberian Snakes fix the design-orbit spin tune Pq to 
^ during the acceleration cycle so that no first-order 
resonances have to be crossed. Crossing first-order 
resonances can lead to a severe reduction of polar- 
ization by an amount described by the Froissart- 
Stora formula (l6| . 

2. Siberian Snakes strongly reduce the influence of en- 
ergy variations on spin motion within a synchrotron 
period 

3. Siberian Snakes reduce the variation of n(z) for par- 
ticles which oscillate vertically and therefore pass 
through horizontal fields which perturb the spin 
motion. 

4. When n changes rapidly during acceleration, the 
adiabatic invariance of Js — n(z) ■ S might be vi- 
olated and polarization would be reduced. It is 
therefore important that Siberian Snakes smoothen 
the changes of ft during the acceleration cycle. 

5. Siberian Snakes can also compensate perturbing ef- 
fects of misaligned optical elements Ha, Ulj but 
the effect of misalignments will not be covered here. 

There is so far no reliable formula for determining the 
number of Siberian Snakes required for an accelerator 
|2lL l22T |. To make things worse, for any given number 
of Siberian Snakes there are very many different possible 
combinations of the snake angles which lead to an energy 
independent closed-orbit spin tune of ^ and to a vertical 
design orbit spin direction no in the accelerator's arcs. 
But so far there has been no reliable formula for deter- 
mining which of these snake schemes leads to the highest 
polarization. 

There used to be a popular opinion that, owing to their 
symmetry, 5 standard choices of the snake angles for 4 
Siberian Snakes are advantageous for HERA-p. These 
choices are not optimal, as will be shown. For reasons 
why these standard schemes were considered useful see 
for example |2^|. RHIC with its two snakes, is operated 
with a similar standard scheme [24| . 

The energy dependence of Pu m in HERA-p produced 
by these 5 snake schemes is shown in Fig. ^ They seem 
to produce rather similar but very low maximum time av- 
erage polarization Pu m in a critical energy regions where 
very strong resonances are excited. The observation of 
such rather small differences in the n-axis for such dif- 
ferent schemes suggests the following detailed investiga- 
tion of the influence of snake schemes. Figure [3 (left) 
shows Pum as computed for linearized spin-orbit motion 
41 for 4 other schemes with 4 Siberian Snakes which were 
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FIG. 1: Pnm of linearized spin-orbit motion for a 2.5<r verti- 
cal amplitude for 4 Siberian Snakes in HERA-p. Each curve 
corresponds to one of the 5 standard choices of snake schemes 
in the following table. The right column shows an obvious 
notation to describe the snake angles in a snake scheme: 



Scheme 


South 


East 


North 


West 


(OfOO) 


0° 


90° 


0° 


0° 


(fOfO) 


45° 


0° 


45° 


0° 


(2f0^0) 


-45° 


0° 


-45° 


0° 


/ 57T 7T 57T 7T \ 

V 8 8 8 S I 


-22.5° 


22.5° 


-22.5° 


22.5° 


( 7T 5-7T 7T 57T \ 
V 8 8 8 8 - 1 


22.5° 


-22.5° 


22.5° 


-22.5° 



found by the so called filtering method 0, |25| , a numer- 
ical search for suitable snake axes. It is apparent that 
large increases in Pu m can result from the choice of a 
suitable snake scheme. 
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FIG. 2: Pnm for 4 different snake schemes which lead to very 
different maximum time average polarization for the standard 
HERA-p lattice. Right: One of these snake schemes which 
were found by numerical optimization of the snake angles. 

HERA-p is the proton ring of an electron-proton col- 
lider and is located above the electron ring. It therefore 
has non-flat sections in which the protons are bend down 
to collide with the electrons. These sections bend spins 
out of the vertical so that no would be non-vertical in 
the arcs. To avoid this, it has been proposed to insert a 
Siberian Snake in each of the 6 vertically bending regions. 
These so called flattening snakes have been assumed in 
all presented computations in addition to the Siberian 
Snakes for which an optimal axis is sought. As a la 
emittance a realistic value of 47rmm mrad was assumed 
throughout this paper. 



III. SPIN-ORBIT-COUPLING INTEGRALS 

For the spin-orbit-coupling integrals in flat rings of 
Eq. n was assumed to point vertically upward. Now 
Siberian Snakes will be included, which rotate all spins 
and also no by tt around some axis in the horizontal 
plane, so that n*o points downward in some sections. This 
changes the sign of 'F Eq. 

It is assumed that there are n Siberian Snakes in the 
ring and that n is even, to make no vertical in the arcs 
of the ring. The azimuth at the position of a snake is 
denoted by 9j and the spin phase advance around the 
vertically upward direction between snake j and j + 1 is 
denoted by ^fj. The spin phase advance after the jth 
Siberian Snake is ^j(0) with Wj(8j) = 0. 

For simplicity 9q — and 9 n +\ = 2-7T is used and the 
spin phase advance from azimuth 8q to the first Siberian 
Snake is 'Jo- 
in the following a Siberian Snake with a snake axis 
which is in the horizontal plane will be referred to as 
a horizontal Siberian Snake and for historical reasons a 
snake which rotates spins around the vertical by some 
rotation angle will be referred to as type III snake. A 
horizontal Siberian Snake with a snake angle tp is equiv- 
alent to a radial Siberian Snake followed by a type III 
snake that introduces an extra spin phase advance of 
a = 2<f H. 

With these notations the spin-orbit-coupling integrals 
for a ring with horizontal Siberian Snakes are 

x J V iM ~ )J( * j(e)+Qj)± * ] d0 • (2) 

In terms of the orbital phase advance $j between snake 
j and j + 1, one obtains 

J=0 

x / ^eH-C-^W+^i^Wldtf . (3) 

A corresponding formula has been used in to in- 
troduce so-called strong spin matching, where Siberian 
Snakes are used to produce a cancelation of spin pertur- 
bations in different FODO cells. 

The spin phase advances between snakes must satisfy 
the condition X)fc=o(~ ) ^fc = to make the closed-orbit 
spin tune independent of energy and the snake angles 
must satisfy the condition Y^=i ( — ) k fk — § to make 
the closed-orbit spin tune equal to \. 
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IV. SNAKE MATCHING IN RINGS WITH 
SUPER-PERIODICITY 

The spin perturbations in different parts of the ring 
can compensate each other when these parts have similar 
spin-orbit-coupling integrals. This is achieved by using 
the Siberian Snakes to adjust the spin phase advances 
in such a way that spin-orbit-coupling integrals of sim- 
ilar parts of a ring cancel each other. In the following, 
the process of finding a snake scheme for which such a 
com pen sation occurs will be referred to as snake match- 
ing [53, . After demonstrating the idea for type III 
snakes, which simply rotate spins around the vertical by 
some fixed angle with little influence on the orbit motion, 
two quite general results will be demonstrated: 

1. A ring with super-periodicity 4 can be completely 
snake matched using 8 Siberian Snakes, i.e. a 
snake scheme can be found for which the spin-orbit- 
coupling integrals are zero due to a complete can- 
celation of spin perturbations in different parts of 
the ring. There are exactly two such possibilities 
which lead to energy independent snake angles. 

2. Such a ring can also be snake matched using 4 
Siberian Snakes. Then, however, the snake axes 
depend on energy and have to be changed during 
the acceleration process. 

3. Energy independent snake matching with 4 
Siberian Snakes can be found when the betatron 
phase advance is appropriately chosen for each of 
the 4 quadrants. 

A. Snake Matching with Type III Snakes for 
Super-Periodicity 4: 

The index y on the spin-orbit-coupling integral and 
on the vertical phase advance and tune will not be in- 
dicated. In any case, the methods for canceling spin- 
orbit-coupling integrals by a special choice of snake angles 
which will now be derived can also be used for transverse 
and longitudinal motion. In this section the notation will 
be further simplified by using the symbols vq and Q to 
denote 2tt times the spin tune and 2ir times the orbital 
tunc. Then for a ring with super-periodicity 4, I ± can 
be computed from 

If = /"^ ! jfc„e i (-* ± *M0 , (4) 
4 Jo 

J± = J±h +e i(-^±Q)/4 + e i2(- I /o±Q)/4 + e i3C-r/o±Q)/4i _ 
4 

Spin transparency requires that I + as well as I~ vanish. 
Thus the bracket in (@J must vanish. This is only possible 
when Q l {- v o±Q)/i is either —1 or i. Choosing the first 
possibility to eliminate / + and the second to eliminate 
I~ , one obtains 

e i(-*o+Q)/4 = _! s e i(-. -Q)/4 = f< (5) 



This leads to the requirement e I( ^/ 2 = i which cannot be 
satisfied in a realistic ring. Therefore, a four-fold repet- 
itive symmetry cannot lead to spin transparency at any 
energy. While the spin disturbance of two quadrants can 
therefore not cancel in I + as well as in I~, one of these 
integrals can cancel whenever the spin phase advance be- 
tween the quadrants is appropriate. 

The situation changes if type III snakes are installed. 
As first found in 25] , type III snakes can improve the spin 
dynamics in HERA-p by increasing Pu m =\{n)\. They 
can be used to manipulate the spin phase advance to 
make the spin-orbit-coupling integrals of different parts 
of the ring cancel. To demonstrate this, 4 type III snakes 
are installed regularly spaced around the ring. 

There are three possibilities for canceling the spin dis- 
turbances between quadrants of the ring. The quadrants 
whose destructive effects cancel are connected by arrows 
in Fig.EI 




FIG. 3: The three possibilities for canceling the depolarizing 
effects of quadrants of a ring with super-periodicity 4. The 
arrows indicate which quadrants cancel 

The spin-orbit-coupling integrals are 

j± = J ±[ 1 + e i(-r/o±Q)/4-Vi ( 6 ) 
4 

_|_ e i2(-r/o±Q)/4-V>i-V2 + e i3(-r/o±Q)/4-i/>i-V>2-V>3] 

where if>j is the spin rotation angle of the type III snake at 
6 = j\. To snake match the ring, I + as well as 1~ must 
vanish. Therefore the bracket on the right hand side has 
to vanish in both cases. A sum of 4 complex numbers 
with unit modulus can only vanish when it consist of two 
pairs of numbers which cancel each other. This is shown 
in Fig. 2| 

The three possibilities of cancelation demonstrated in 
Fig. are given by the following three sets of equations: 

1. (-I/Q ± Q)/4 - fa = 7T 

and (— 1/ ± Q)/4 — ip 3 = n 

2. 2(-i*±Q)/4-Vi-V>2 = 7T (?) 
and ip3 = ipi 

3. s{-v ±Q)/4-i> 1 -il> 2 -il> a = 7T 
and (— v a ± Q)/4 — ^2 = 7r 

The symbol = indicates equivalence modulo 2ir. To 
snake match, one of these three conditions has to hold 
for (—vq + Q), which lets / + vanish and another of the 
conditions has to hold for (—^0 — Q), which lets /~ van- 
ish. / + and / _ cannot vanish due to the same condition 
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FIG. 4: Four complex numbers with modulus one can only 
add up to zero when they consist of two pairs which individ- 
ually add up to zero 

if restrictions on the allowed orbital phase advance Q are 
to be avoided. There are therefore three possibilities: 

1. J + = due to condition 2 and I~ = due to 
condition 3 requires 

il>2 = 7r + 2(-^o + Q)/4-Vi , ^3 = ^1, (8) 
il>2 = ir+(-vo-Q)/4, ifo=2(-va-Q)/4-tl> 1 . 

The first and the third of these equations require 
that tpi = (— ^o+3Q)/4 whereas the second and the 
fourth equations require that ipi = (—i>o — Q)/4. 
These two requirements are in general not compat- 
ible and the ring cannot be made spin transparent 
in this way. 

2. I + = due to condition 1 and I~ = due to 
condition 3 requires 

1>l = * + {-vo + Q)/4 , ip 3 = ih. , (9) 
02 = n+(-is Q -Q)/4, V'3 = 2{-vo - Q)/4 - Vi • 

The first and the last of these equations together 
require 03 = tt— (3Q— Vq)/A. This is in conflict with 
the second equation. Thus this way also cannot 
lead to a spin transparent ring. 

3. I + = due to condition 1 and I~ = due to 
condition 2 requires 

i>i = ^+(-^o + Q)/4, ^3 = 01, (10) 

02 = 7T+2(-I/Q -Q)/4- V>1 i ip3=ipi- 

These 4 equations are compatible and lead to ipi = 
^3 = 7T + (-1* + Q)/4 and V2 = (-^b - 3Q)/4. 

The type III snake at I = has the rotation angle 
i/>4 which is chosen in such a way that the closed-orbit 
spin tune of the ring does not change due to the snakes, 
i.e. ipi + -02 + V'3 + "04 = 0- The required rotation angles 
are then 

■01 = 03 = 7T H , (11) 

, o i/ + 3Q , o Q - 3 - fo 

02 = — , 04 = i ■ 



Obviously a change in sign of Q leads to I + — due 
to condition 2 and to I - = due to condition 1. There 
are therefore exactly two possibilities for making a ring 
with super-periodicity 4 spin transparent by means of 4 
type III snakes. These possibilities are shown in Fig. [S] 
However, the scheme of 4 type III snakes presented here 




FIG. 5: The only two ways to snake match a ring with super- 
periodicity 4 by 4 type III snakes. The number [x] denotes 
x mod 2n and lies in [0, 2ft). The vertical tune times 2ir is 
denoted by Q and vo = G^2n 

cannot be a practical snake scheme, since it does not 
make the closed-orbit spin tune independent of energy. 
But it illustrates how type III snakes can be used at fixed 
energy to make spin perturbations from different parts 
of the ring cancel each other. This feature can then be 
used in combination with the Siberian Snakes which are 
installed to make the closed-orbit spin tune independent 
of energy. 

B. Snake Matching with Type III Snakes for 
Super-Periodicity 4 and Mirror Symmetry: 

In particle optical systems, mirror symmetries are of- 
ten used to cancel perturbative effects [2|| |3(], 0S1 H3 • 
Therefore it is interesting to see whether mirror symme- 
try can lead to vanishing of spin-orbit-coupling integrals 
when 4 Siberian Snakes are installed at the symmetry 
points of the ring. If one super-period is mirror symmet- 
ric, then 

If = r 4 A v e I( "* ± * )dfl . ( 12 ) 
8 Jo 

If = If + (/±)* c i(— o±Q)/4 _ (13) 

4 8 8 

With E ± = (-u ± Q)/4 one obtains for the complete 
ring 

7± = (If + (If )*e is± ) • (1 + c iS± + e i2S± + e i3S± ).(14) 

8 8 

Type III snakes will not be considered further, since 
a horizontal Siberian Snake with snake angle tp can be 
decomposed into a radial Siberian Snake and a type III 
snake with rotation angle 2<p. 
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C. Snake Matching of Siberian Snakes with Fixed 
Axes at All Energies for Super-Periodicity 4: 

Thus snake matching the ring with super-periodicity 4 
is not influenced by the fact that the ring might have a 
mirror symmetry since the bracket in (|14f> is equivalent to 
the corresponding bracket in Q for rings without mirror 
symmetry. 

a. Schemes with 4 Snakes: For 4 horizontal Siberian 
Snakes the spin-orbit-coupling integral in Q is 



7± = 



+e i(-* ±* ) I fcye i(*i+"i±*i) d fl 



(15) 



_|_ e i(-*o+ai+*i±(*o+*i)) 

' fc„e i( -* 2 - Q2± * 2) d0 



_|_ e i(-*o+ai+*i-«2-*2±( < ! , o+*i+*2)) 



k i(*3+a 3 ±*3) d6 , 



For a ring with super-periodicity 4 and with 4 equally 
spaced horizontal Siberian Snakes one obtains with i>q = 
*(27r), = *(2tt)/4, and ^ = Q/4 the relation 



J? 



(16) 



J± _ jf(l _|_ c i(-*o+cn+*l-c«2±( < t , o+*l))^ ^7) 

4 

+(If )*e i (~* 0+Ql± *°) 

4 

x (1 + e i (* 1_Q2 ~* 2+ct3± (' I ' 1+ * 2 ^) 

_ ji/j + e i(ai-a 2 ±2Q/4)^ 



j((-v ±Q)/i+ ai ) 



(1 + gi(-a2+a 3 ±2Q/4h 



Spin transparency of the ring is therefore obtained when 

cti - a 2 ± 2Q/4 = 7T and - a 2 + a 3 ± 2Q/4 = tt . (18) 

This cannot be achieved in general since the conditions 
cyi — a-j + 2Q/4 = 7r and a.\ — a 2 — 2Q/4 = n have to be 
satisfied simultaneously, which implies Q = 0. 

In the case of mirror symmetry in the ring if and IT 

4 4 

are related by ifH3|) . 

if = If + (/±)*e i (-"o±Q)/4 ; ( 19) 

4 8 8 

(If)* = {If )* + If c-^-^TQ)/* 

4 8 8 

= If e-K-"oTQ)/4 . (20) 

4 

With mirror symmetric quadrants the spin-orbit- 
coupling integral then simplifies to 

j± _ j±Q + e i(ai-a 2 ±2Q/4)\ 

4 

+ 7f (1 + e i (-«2+« 3 ±2Q/4)\ e i(a !l ±2Q/4) (21) 



and again this additional symmetry does not simplify the 
compensation of the spin-orbit integrals. 

b. Schemes with 8 Snakes: The same procedure can 
now be repeated with 8 snakes. For that purpose 4 more 
horizontal Siberian Snakes are placed at the locations 
j7r/2 + A9, j e {0, 1, 2, 3}. In terms of the integrals 



AH 



It = / kyj^^M, 

JA6 

the spin orbit coupling integrals of J2J are 



(22) 
(23) 



j± _ /±(1 -|- e i(-*o+ai+*i-c(2±(*o+*i)) _|_ (24) 

_|_ c i(-*o+ai+*i-a2=F----"ii-2±(*o + ...+*, 1 -3))^ 

_|_ j± e i(-* +ai±*o)( 1 _|_ e i(*i-a 2 -*2+a3±(*i+*2)) 

+ ... 

_|_ c i(*l-Q2-*2+Q3±--- + Ct„-l±(*l + --- + * n -2))^ 

If there is an additional mirror symmetry and the snakes 
are all placed in the symmetry points, (|13fl implies lf~ = 
(Iq )*e'( -!yo±( 2)/ 4 , which again does not lead to simpli- 
fications. The complete spin phase advance of the ring 
is X)j=o( — y i^j + a i) = ""• Since this phase advance 
is required to be independent of energy, X)J=o(~ 
has to vanish. Because of the super-periodicity this re- 
quires U'o = 4*1, and all the spin phases 4* j in the equa- 
tions (|25|l cancel. Then in terms of the difference angles 
Ajk = ctj — ctk, spin matching the ring requires 

l + c i(±Q/4+A 12 ) + e i(±2Q/4+A 12 + A 34 ) 

+ e i(±3Q/4+A 12 +A 3 4+A 5 6) _ j (25) 

1 + e i(±Q/4-A 23 ) +e i(±2Q/4-A 2 3-A 45 ) 

+ e i(±3Q/4-A 23 -A 4 5-A 6 7) _ q _ (26) 



Sets of 4 complex numbers with modulus 1 can only add 
up to zero by the three schemes shown in Fig. The 
equations (|25ll and H26|l have the same structure as the 
matching conditions of equations JfjJ) and the relations 
(|12fl can therefore be used to obtain the following two 
ways to satisfy (I25|l : 

A 12 = A 56 = tt - Q/4 , A 34 = 3Q/4, (27) 
A12 = A 56 = tt + Q/4 , A 34 = -3Q/4. (28) 

Equation l|28|) resulted from reversing the sign of Q in 
1)27(1. There are also exactly two possibilities for solving 

(H3>, 

A 23 = A 67 = tt + Q/4 , A 45 --3Q/4, (29) 
A 23 = A 67 = tt - Q/4 , A 45 = 3Q/4. (30) 

There are now 4 possibilities to snake match the ring; 
these are obtained by combining the equations l|27|l fc (|29|l . 
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O&lUnj, or pS jl fcP ^I . where the last two 

possibilities also result from the first two by reversing the 
sign of Q. FigureEJshows how parts of the ring cancel the 
depolarizing effects of other parts in these snake matching 
schemes. 




FIG. 6: The two possibilities by which the individual parts 
of the ring cancel the spin-orbit-coupling integrals. Changing 
the sign of Q leads to two corresponding snake schemes. The 
symbols Iq , Iq , and 1^ indicate which part of the spin- 
orbit-coupling integrals are canceled 

Since only differences in the snake angles appear, one 
of the angles can be chosen arbitrarily. This then fixes all 
other snake angles. Here ax = is chosen for simplicity. 
This leads to the following possibilities: 

Combination of the equations (|27[) and l|29|) : 

Qti = , Qf2 = TT + Q/4 , Qf3 = , (31) 

a 4 = -3Q/4, a 5 = 0, a 6 = tt + Q/4 , (32) 
a 7 = 0, a 8 = 7r + Q/4. (33) 

Combination of the equations l|27(l and (|30|l : 

ax = 0, a 2 = 7r + Q/4, a 3 = 2Q/4 , (34) 
a 4 = -Q/4, a 5 = -4Q/4, a 6 = tt - 3Q/4 , (35) 
a 7 = -2Q/4, a 8 =7r-Q/4. (36) 

The values for as were obtained from the requirement 
ag — a 7 + a§ — «5 + a 4 — a 3 + a 2 — a 4 = 7r. The last snake 
scheme can be simplified by decreasing all snake angles 
by 2Q/4, leading to 

a x = Q/4, a 2 = 7r + 2Q/4, a 3 = 3Q/4 , (37) 
a 4 = 0, a 5 = -3Q/4, a 6 = tt - 2Q/4 , (38) 
a-j = —Q/4, as = 7r . (39) 

These snake schemes are shown in Fig. where account 
has been taken of the fact that the actual angle between 
the snake's rotation axis and the radial direction is a/2. 
Furthermore advantage has been taken of the fact that 
the angle a/2 only needs to be known modulo tt. 

Here it is very important to note that the snake angles 
are independent of v = G7 and therefore that a snake 




1= n 
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FIG. 7: Two of the 4 ways to snake match a ring with super- 
periodicity 4 using 8 horizontal Siberian Snakes. The number 
q denotes Q/8 mod tt and lies in [0, 7r). The vertical tune 
times 2tt is denoted by Q and v — G^2tt. The other two 
possible snake schemes are obtained by reversing the sign of 
Q. When all snake angles are increased by the same amount, 
then the ring remains spin transparent. Note that the snake 
angle is independent of v and thus of energy 



match has been achieved for all energies. With 4 Siberian 
Snakes such an energy independent snake match is not 
possible in a four-fold symmetric ring. 

One could repeat the same procedure for a layout with 
6 horizontal Siberian Snakes or with combinations of, for 
example, 6 horizontal Siberian Snakes and two type III 
snakes. Due to their three-fold symmetry, this could be 
of special interest for the RHIC rings. 



V. SNAKE MATCHING HERA-P 
A. Schemes with 4 snakes: 

When the spin-orbit-coupling integrals starting at an 
azimuth #0 are minimized, the opening angle of the in- 
variant spin field at #0 for the approximation of linear 
spin-orbit motion is also minimized. In fact, Pu m = 1 
in linear approximation if / = since then the spin 
motion is decoupled from the orbit to first order. 

It has been seen from the example of a ring with super- 
periodicity 4 that 8 Siberian Snakes can be used to snake 
match the spin-orbit-coupling integrals to zero at one az- 
imuth of the ring for all energies. This could be achieved 
since the snake angles were used to adjust the spin phase 
advances in such a way that perturbations in one part of 
the ring were compensated by identical perturbations in 
one of the identical super-periods of the ring. 

HERA-p has 4 identical arcs separated by 4 not identi- 
cal straight sections. Due to the lack of four-fold symme- 
try it is in general not possible to find snake angles which 
completely compensate all spin-orbit-coupling integrals. 
However, the 4 identical arc sections of HERA domi- 
nate the spin-orbit resonance strength of vertical motion. 
Thus it would make sense to arrange that the perturbing 
effect of these arcs cancel each other. A first step in this 
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direction is a symmetrization of the quadrants by mak- 
ing the spin phase advance in all of the straight sections 
identical. This can be done by the insertion of two more 
flattening snakes. 

The spin-orbit-coupling integrals from the first regular 
FODO cell to the last FODO cell of a regular arc in 
HERA-p will be denoted by and I~ and the azimuths 
of the beginnings of the 4 regular arcs as 9\, 82, #3, and 84. 
The central points of the South, West, North, and East 
straight sections are denoted by S, W, N, and E. The 
spin phase advances between the arcs are compensated 
using, the snake angles (fE, fN, and ipw- The spin phase 
advance between 6j and 6j is denoted by These 
notations are indicated in Fig. |S] 




FIG. 8: The spin phase advance from the beginning of one 
regular arc to the beginning of the regular arc on the opposite 
side of the ring. Large arrows indicate Siberian Snakes, mall 
arrows indicate flattening snakes. 

With Siberian Snakes in each of the straight sections, 
the spin phase advance from 6\ to 83 is given by ^13 = 
"Jib — 2(fiE — ^en + 2<Pn + ^N3- The spin phase advance 
is identical in all quadrants of the ring. The total spin 
phase advance is then solely determined by the snake 
angles and is therefore independent of energy: ^13 = 
2(ifiN — <Pe) an d ^24 = 2(<y9vF — <Pn)- The orbital phase 
advance &(8 3 ) — <f>(6>i) also does not depend on energy. 
For simplicity, $(0?) — $(#i) will now be denoted by $^ . 

The spin-orbit-coupling integrals at the South interac- 
tion point then contain the following contributions from 
the 4 regular arcs: 

/ a + cs = i+c i (-* si +*si)(i + e i[2(^- v »)+* 1 3] ) ( 40 ) 

+ (/ _ )*e i ' 2¥ ' E_ * S - B + * E2 + * S2 ' (1 + e i [ 2 (V»'-VN) + *24]) 

J- = /- e i(-*si-*si)(i _|_ e i[2(vE-pjv)-*i3]) (41) 

I ( f+ N l* p i(2ipE-*SE+'I'E2-*S2)('-| I p i[2(>PW->PN)-$2i]\ 



This shows that it is always possible to cancel one of 
the spin-orbit coupling integrals by choosing the snake 
angles so that the spin perturbation produced in one of 
the arcs is canceled by the arc on the opposite side of 
the ring. Since |/ + | and |/ _ | are different, neighboring 
arcs can in general not compensate each other when 4 
Siberian Snakes are used. 

It is however possible to use the eight-snake scheme 
found for symmetric lattices. The two special 8 Siberian 
Snake schemes which lead to an energy independent 
snake match in a ring with super-periodicity 4 will not 
spin-match HERA-p completely, but the spin perturba- 
tion from the arcs, which are the dominant perturbation, 
will be compensated exactly. This possibility of having 
a set of Siberian Snake angles which do not have to be 
changed with energy and which lead to a tightly bundled 
invariant spin field is on the one hand very attractive; on 
the other hand it requires 8 Siberian Snakes of which 4 
would have to be installed at the centers of the HERA-p 
arcs, where technical requirements of moving cryogenic 
feed-throughs and super-conducting magnets would be 
very costly. If possible, a four-snake scheme should there- 
fore be found. 

Whereas it was shown below (|18|) that a four-snake 
scheme cannot cancel both spin-orbit-coupling integrals 
in a ring with super-periodicity, a corresponding cancela- 
tion of the spin perturbation due to the arcs in HERA-p 
can nevertheless be achieved since the orbital phase ad- 
vances between the arcs can be manipulated individually, 
while these four phase advances are equal for a lattice 
with super-periodicity 4. 

To cancel both spin-orbit integrals in 1421) . 4 phase 
factors have to be — 1. This requires 

2(<Ae - f N ) + $13 = 7T , (42) 
2(ip E - <Pn) - $13 = T , (43) 

2((fl W - ip N ) + $24 = 7T , (44) 
2(if W - ip N ) - $24 = 7T ■ (45) 

Subtraction of the first two equations leads to the re- 
quirement that the betatron phase advance from 8\ half 
way around the ring to #3 is an odd or even multiple of 
7T. The same is true for the phase advance from 82 to 
84. Correspondingly, the spin phase advance over these 
regions has to be an odd multiple of 7r when the orbit 
phase advance is an even multiple and vice versa. With 
a rather benign change of the vertical optics in HERA-p 
which does not change the vertical tune, the contribution 
of the regular arcs to both spin-orbit-coupling integrals 
can be canceled, even in a four-snake scheme. 

The snake scheme (Of §f )o has $13 = and ^24 = 0. 
For this snake scheme, the betatron phase advances from 
6*i to #3 and from 82 to 84 were adjusted to be odd mul- 
tiples of 7r. This change of the optics is indicated by the 
index o in the notation of the snake scheme. The maxi- 
mum time average polarization Pu m is plotted (blue) in 
Fig. El for the complete range of HERA-p momenta (top) 
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and for the critical momentum regions above 800 GeV/c that the snake match is independent of energy since then 
(bottom). As a comparison, Pu rn for a standard snake the nth octant's contribution JfLj to the spin-orbit cou- 



scheme (^0^0) (red) is also shown. 

The complete snake match of the arcs in HERA-p in 
deed eliminates all strong reductions of Pu m over the 
complete momentum range. Nonlinear effects will be an- 
alyzed later, but as far as the linear effects are concerned, 
this snake matched lattice of HERA-p would be a rather 
promising choice for the acceleration of polarized proton 
beams. 




310 820 830 

p (GeV/c) 



FIG. 9: Improvement of Pu m by matching 4 snake angles and 
the orbital phases. The snake arrangement is (0^^--|)o (blue). 
As a comparison Pu m from linearized spin-orbit motion is 
shown for the same HERA-p optics with a (^0j0)o snake 
scheme (red) 



B. Schemes with 8 snakes: 

Although, for the reasons explained, eight snakes are 
not very practical for HERA-p, significant improvements 
would be possible if 8 snakes were used, as will now be 
shown. The snake matching schemes of Fig. {7\ are suit- 
able for rings with super-periodicity 4. They are there- 
fore not directly applicable to HERA-p. However, the 
cancelation schemes of Fig. [H] can still provide a guide to 
snake matching HERA-p with 8 snakes. 

The 8 Siberian Snakes are placed in the straight sec- 
tions and into the centers of the arcs, so that the hori- 
zontal angle between adjacent Siberian Snakes is 45° and 
therefore tyj = for all j £ {0 ... 7}. This ensures 



pling integral 7 can be compensated by the 2nd neigh- 



bor's contributions 7„ 
the phase factor in 



In particular the \Pj cancel in 



(46) 



,![(-!)"(*„_!- 



-*„+a„ + i)±(*„-i+*„)] 



Now a snake scheme for HERA-p is sought that is 
guided by the cancelation scheme of Fig. |SJ (left). When 
the first octant's contribution 7q to a spin-orbit coupling 
integral 7 ± is to be compensated by that of its 2nd neigh- 
bor the phase factor in 



a i[-*o+ai + , I'i 



-a 2 ±(*o + *l)] t± 
1 2 



(47) 



should be 
factor in 



-1. The same should be true for the phase 



-*4+a 5 +*5-Q6±(*4 



M&<s)]j± 



(48) 



These conditions are satisfied for the superscript + as 
well as — when all the snake angles ctj are zero for 
j € {I, 2, 5, 6} and the betatron phase advances are cho- 
sen appropriately. Similar conditions arise for the phase 
factors involved in matching 7i against 73 and in match- 
ing 75 against I?. They are also —I when the betatron 
phases are chosen appropriately and when 0:3 and aj are 
zero. Since ao and 014 do not appear in these matching 
conditions, they can be chosen freely. But for a design 
orbit spin tune of 0.5, one has to choose ao + an = tt. 
Here CV4 — has been chosen and this snake scheme is 
then characterized as 0000000) o. 

In the case of HERA-p spin-orbit integrals of complete 
octants cannot cancel, but the contribution of the regu- 
lar arc sections can cancel. Therefore the betatron phase 
advance is not chosen to be $0 + ^1 = tt between Siberian 
Snakes, but rather the betatron phase advances between 
the beginning of the regular arcs of the first and the sec- 
ond quadrant are chosen in that way, i.e. $12 = 7r with 
the notation of Fig. |H1 In this configuration the regular 
arc of the first octant cancels that of the third, i.e. the 
regular arc's contribution to 7 2 cancels that to 7 . But 
also the regular arc of the second octant cancels that of 
the fourth since the betatron phase advance between the 
centers of the first and the second arc is then also tt mod 
2-7T, i.e. the regular arc part of 73 cancels that of Ii . Simi- 
larly the octants of the third quadrant cancel those of the 
fourth quadrant. In the snake scheme (^0000000)0 the 
phase advance over the East straight section was changed 
by 2tt x 0.1028 to have $12 = 2tt x 8.5. The phase ad- 
vance over the West was changed by 2ir x 0.0208 to have 
$34 = 2tt x 7.5. In order to have the same betatron phase 
advance in the North and the South straight sections and 
to keep the total tune constant, the linear optics of the 
North and of the South straight sections were modified. 

This cancelation scheme does not agree completely 
with that of Fig. ©(left) since in (§0000000)o the super- 
scripts + and — have been dealt with simultaneously so 
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that Iq no longer has to compensate If . Since $12 = 7r 
and $34 = 7r, it is appropriate to take q = tt/2 for the 
phase advance of one octant of the scheme in Fig. (left), 
for which one then finds a close similarity to the snake 
scheme (^0000000)o. The resemblance would be even 
closer if the phase advance over the North straight sec- 
tion had been changed so that $23 = and if = n, 
ao = had been chosen. This would require an addi- 
tional change of $41 to adjust the orbital tunc. 

A snake scheme for HERA-p which resembles 
Fig. (right) can also be found. In this cancelation 
scheme the phase factors in the following sums have to 
be -1: 



(49) 
(50) 
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FIG. 10: Improvement of Pu m by matching 8 snake angles and 
the orbital phases. The snake arrangement is (^0000000)o 
[blue) and (^abc0-c-b-a)o (green). As a comparison Pn m from 
linearized spin-orbit motion is shown for the snake and phase 
advance matched four-snake scheme (0-|-|-|)o (red) 

Again complete octants cannot cancel, but to cancel 
the contributions of regular arc sections, the beginning 
of a regular arc is chosen as the starting point for the 
compensation. Now the <E>„ for all even n describe the 
betatron phase advance from the beginning to the center 
of the regular arc, and they are therefore equal. For the 



choice 0:4 = all phase factors are —1 if and only if 
$ 5 = <fr 1 and 

Oil = $23 j "2 = $13 , «3 = $14 1 (57) 

a 5 = -a 3 , a 6 = -a 2 , a 7 = -ai . (58) 

This snake scheme is referred to as (^abc0-c-b-a)o. The 
condition $5 = $1 was satisfied by changing the betatron 
phase advance over the East and West straight section 
without changing the vertical tune. 

Linearized spin orbit motion leads to a very favorable 
Pu m for both eight-snake schemes as shown in Fig. 
where it is compared to the Pu m of the snake and phase 
advance matched four-snake scheme (05-|-?)o. 



VI. HIGHER-ORDER RESONANCES AND 
SNAKE SCHEMES 

At the critical energies, where the maximum time av- 
erage polarization is low during the acceleration process, 
linearized spin-orbit motion does not describe spin dy- 
namics well. The spin motion is influenced by several 
overlapping resonances in these regions and the single res- 
onance approximations jl6l |33| can also not be applied. 
Thus the simulation results obtained with these computa- 
tionally quick techniques should always be checked with 
more time consuming non-perturbative methods if pos- 
sible. This is also true for the snake-matched lattices of 
HERA-p with large Pu m , even though they avoids large 
variations of the invariant spin field n(z) over the phase 
space of the beam in linearized spin-orbit motion. When 
first-order effects are canceled, the higher-order effects 
become dominant and the quality of the snake-matched 
lattice of HERA-p can only be evaluated with higher- 
order theories. 

Until 1996, when stroboscopic averaging 01 was in- 
troduced, there was no non-perturbative method of com- 
puting the n-axis at high energy in proton storage rings, 
where perturbative methods are usually not sufficient 
[l4| . In addition, the method of anti-damping was de- 
rived [34|], which also computes n(z) non-perturbatively 
and which can be faster when the n-axis is required for a 
range of phase space amplitudes. Both methods of com- 
puting the invariant spin field are implemented in the 
spin-orbit dynamics code SPRINT, by which also the 
amplitude-dependent spin tune v(J) can be computed 
once n(z) is known. Since stroboscopic averaging and 
anti-damping are based on multi-turn tracking data, they 
are applicable to all kinds of circular accelerators and 
they are especially efficient for small rings and for simple 
model accelerators. Nonlinear motion for orbit and spin 
coordinates can esily be included [35ll36| . 

Subsequently another non-perturbative algorithm for 
computing n(z) and v(J) has been derived (3?| ■ It is 
called SODOM-2 since it was inspired by the earlier al- 
gorithm SODOM ^3 which for convergence required the 
angle between n and no to be small. With some routines 
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provided by K. Yokoya, SODOM-2 was incorporated into 
the program SPRINT [ll HI and leads to results 
which agree very well with those of stroboscopic aver- 
aging. For motion in one degree of freedom, SODOM-2 
is often faster than stroboscopic averaging, especially for 
large rings like HERA-p where particle tracking is rela- 
tively time consuming. But for orbit motion in more than 
one degree of freedom or in the vicinity of spin-orbit res- 
onances, SODOM-2 becomes exceedingly slow and then 
stroboscopic averaging and anti-damping are needed. 

To check whether the improvements of spin motion ob- 
tained in the framework of linearized spin-orbit motion 
survive when higher-order effects are considered, Pu m 
and v has been calculated. The result for one of the 
standard Siberian Snake schemes which used to be con- 
sidered advantageous by popular opinion is shown for 
the South interaction point of HERA-p in Fig. ^2 R 
has four Siberian Snakes in the (?0f 0) scheme. Some 
of the features of Pu m were already revealed by lin- 
earized spin-orbit motion in Fig. ^ Now many higher- 
order resonances are revealed, causing strong reduction of 
Piim and there are corresponding strong variations of the 
amplitude-dependent spin tune v |4pL l4l| . Strong reso- 
nances occur especially in the critical energy region where 
linearized spin orbit motion in Sect. ^already indicated a 
very small Pu m due to a coherent spin perturbation in all 
regular FODO cells. Many higher-order resonances over- 
lap in these critical energy regions of Fig. ^] (top) where 
large spin tune jumps can be observed in Fig. 1111 (bot- 
tom). The strongest spin tune jumps occur in the criti- 
cal energy regions, mostly at the second order resonance 
v = 2Q y which is indicated by the top line |4^. 




p(GeVfc) p(GeV/c) 

FIG. 11: Pum and v for particles with a vertical amplitude 
corresponding to the 2.5a emittance in the HERA-p lattice of 
the year 2002 with the (^0^0) scheme. Top: The maximum 
time average polarization Pn m for the complete acceleration 
range (left) and for the critical energy range above 800 GeV/c 
which has to be crossed when accelerating to the proposed 
storage energy of 870 GeV/c. Bottom: the corresponding 
amplitude dependent spin tune v(J y )- The second-order res- 
onances v = 2Q y and v = 1 — 2Q y are indicated (red) 

Figure 1121 also shows Pu m and v for higher-order 
spin dynamics in the snake-matched and phase-advance- 
matched HERA-p ring with 4 Siberian Snakes. While 



the overall behavior of Pu m over the complete accelera- 
tion range of HERA-p looks similar to the result obtained 
with linearized spin-orbit motion, which was displayed in 
Fig. El higher-order effects become very strong at high 
energies, especially in the vicinity of the critical ener- 
gies where perturbations of spin motion in each FODO 
cell accumulate. The spin tune spread at momenta be- 
low 400 GeV/c is small and higher-order effects seem to 
be benign even at these critical energies. A comparison 
of Figs. ITT| and 1121 hello shows that the special scheme 
obtained by matching orbital phases and snake angles 
would be a very good choice for HERA-p up to 300 or 
400 GeV/c. For linearized spin-orbit motion, this snake- 
matched scheme does not produce a strong reduction of 
Piim at any critical energy, but the higher-order effects 
become very pronounced at some top energies of HERA- 
p. Nevertheless, the advantage over other snake schemes 
becomes clear in figure Fig. ^] (bottom) where the am- 
plitude dependent spin tune v is shown. It comes close 
to a second order resonance at fewer places and does not 
exhibit spin tune jumps which are as strong as those in 
previous figures of v. 

In this snake-matched scheme, the influence of higher- 
order effects can be seen very clearly, because the first- 
order effects have been matched to be very small. The 
analysis of v shows that completely snake matching the 
spin perturbations in the arcs of HERA-p with 4 Siberian 
Snakes is advantageous, even though dips of Pu m due to 
higher-order resonances can be observed at high energies. 
Even around 300 GeV/c there are resonant dips of Pu m 
in Fig. E| (middle) but they are less pronounced than 
those in Fig.^Jso that the snake matched scheme should 
be very advantageous in the complete energy range of 
HERA-p. 



VII. POLARIZATION REDUCTION DURING 
ACCELERATION 

It should be noted that the destructive spin tune jumps 
at second order resonances disappear completely when 
HERA-p is simulated without its non-flat regions. This 
is due to the fact that a large class of resonances are not 
excited at all in mid-plane symmetric rings 0, l43l | . To 
reduce these perturbations, the East region of HERA-p 
will now be simulated as flat since the HERMES exper- 
iment located in this region does not require that the 
proton beam is on the level of the electron beam. 

When a particle is accelerated across the critical mo- 
mentum region from 800 to 806 GeV/c with a typical 
acceleration rate of 50 keV per turn, the adiabatic in- 
variance of Js — n ■ S can be violated and the level of 
violation will depend on the orbital amplitude and the 
snake scheme. This violation is illustrated in the graphs 
in Fig. 1131 which, for three different snake schemes, show 
the average spin action Js at 806 GeV/c which had ini- 
tially Js = 1 at 800 GeV/c before acceleration. 

The change of Js in the critical energy region depends 
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FIG. 12: Pum and v for a 2.5<r vertical amplitude after the 
betatron phase advance between opposite regular arc struc- 
tures was adjusted to be an odd multiple of tv in the (0? S-^o 
scheme. The resonances v — 2Q y and v = 1 — 2Q y are indi- 
cated (red) 



on the initial phase space angle so that if Js had been 
computed only for one particle, it could by chance have 
had an angle variable for which Js does not change al- 
though it would have changed for other points with the 
same vertical phase space amplitude. To avoid such a 
chance effect which gives the impression that Js is in- 
variant, three particles were accelerated and the average 
Js is displayed in Fig. ^] 

At small phase space amplitudes, Js is nearly invari- 
ant and therefore J$ = 1. For each of the three snake 
schemes, there is a phase space amplitude J ym ax above 
which Js < 1 and the regions of the beam with an am- 
plitude above J ym ax lead to a reduction of the beam's 
polarization during the acceleration process. 

For the standard snake scheme (f 0^0), only the part 
of the beam with less than Itt mm mrad vertical am- 
plitude can remain polarized. For the filtered scheme 
(tTTt) 9 ' phase space amplitudes up to 4tt mm mrad 
are allowed. Finally the snake matched scheme (jOjj)o 
gives the most stable spin motion and Fig. ll3l shows that 



vertical amplitudes of up to 87r mm mrad are allowed. 
This shows that the snake matched scheme is superior 
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FIG. 13: Average spin action Js at 806 GeV for particles 
starting with Js = 1 at 800 GeV/c for three different snake 
schemes. Left: (fOfO), Middle: (^||f)o, Right: 
(0^??)o scheme. Particles with an amplitude above 1 (left), 
4 (middle), and 8 (right) lead to a reduction of polarization 
when the beam is accelerated through this critical energy re- 
gion 

to the other four-snake schemes studied here. It stabi- 
lizes spin motion for 10 times larger phase space am- 
plitudes than some other snake schemes. Nevertheless, 
Stt mm mrad is not enough to allow high polarization at 
top energies for today's emittances in HERA-p. 

Matching 8 Siberian Snakes in HERA-p lead to two 
snake schemes with very high Pu m and very small spin 
tune spread. But to demonstrate that it is possible to 
further stabilize spin motion in HERA-p by such schemes, 
Fig. O shows the vertical phase space amplitudes for 
which Js remains invariant. The more effective of the 
two snake scheme stabilizes spin motion up to a vertical 
amplitude of 147r mm mrad. 
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FIG. 14: The average spin action Js at 806 GeV for particles 
which started with Js = 1 at 800 GeV/c for the 2 different 
snake matched eight-snake schemes (-|0000000)o (left) and 
(^abcl-c-b-a)o (right). Particles with an amplitude above 13 
(left) and 15 (right) lead to a reduction of polarization when 
the beam is accelerated through this critical energy region 

These results for the various snake schemes are col- 
lected in Fig. El where it becomes clear that snake 
matching with 4 and especially with 8 snakes leads to 
a significant improvement. In HERA-p it does not suf- 
fice to avoid a reduction of Js for particles with less than 
147T mm mrad amplitude. It would therefore be very 
helpful to use electron cooling in PETRA 0, 14f3 ? ] so 
as to reduce the emittance in HERA-p and to allow for 
an acceleration without loss of polarization for most par- 
ticles in the beam. 
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FIG. 15: The average spin action Js at 806 GeV for par- 
ticles which started initially with Js = 1 at 800 GeV/c. 
Violet: The standard scheme which stabilizes spin motion 
for particles within 1-k mm mrad, Cyan: the filtered four- 
snake scheme stabilizes within 4n mm mrad, Red: the snake- 
matched four-snake scheme stabilizes within 8-7T mm mrad, 
Green: the snake-matched eight-snake scheme which stabi- 
lizes within 137T mm mrad, Blue: and the snake-matched 
eight-snake scheme which stabilizes within 147r mm mrad of 
vertical phase space amplitude 



The excellent performance of the two schemes with 8 
Siberian Snakes is due to much smaller oscillation of the 
amplitude dependent spin tune during the acceleration 
process, and the destructive second order spin-orbit res- 
onances indicated in Fig. ^3 (bottom and top-right) are 
hardly encountered when these snake schemes are cho- 
sen. Correspondingly Pu rn only drops to small values at 
very few energies in Figs. ITfjl (center and top- left). 

Contemplating all these results, it can be concluded 
that it is not possible to give a simple formula for the 
number of snakes which are required for a given acceler- 
ator since different snake schemes with the same number 
of snakes lead to very different stability of spin motion. 
It has even been shown that 8 snakes are not necess arily 
better than 4 snakes for the non-flat HERA-p ring |47| . 
In the end detailed evaluation is needed. 



Acknowledgments 

Desmond Barber's careful reading and improving of the 



K. Steffen: Strong spin matching with and without 
snakes, a scheme for preserving polarization in large ring 
accelerators. DESY-88-068 (1988) [12 
E. D. Courant, S. Y. Lee, and S. Tepikin. AIP Conference 
Proceedings 145, page 174, New York, 1985. American 
Institute of Physics. 

D. P. Barber, K. Heinemann, G.H. Hoffstaetter, and [13 
M. Vogt: The permissible equilibrium polarization dis- 
tribution in a stored proton beam. In ICFA Proceed- 
ings Quantum Aspects of Beam Physics, World Scientific 
(1998) [14 
G.H. Hoffstaetter: Aspects of the Invariant Spin Field 
for High Energy Polarized Proton Beams. Habilitation, 
Darmstadt University of Technology, (January 2000) [15 
M. Vogt: Bounds on the maximum attainable equilib- 
rium spin polarization of protons in HERA. Dissertation, 
Universitdt Hamburg, DESY-THESIS-2000-054 (Decern- [16 
ber 2000) 

M. Vogt, D.P. Barber, and G.H. Hoffstaetter: Introduc- 
tion to the invariant spin field and update on simulations [17 
with polarized protons at HERA-p. In DESY-HERA-98- 
04 (1998) 

G.H. Hoffstaetter, H.S. Dumas, and J.A. Ellison: Adia- [18 
batic invariants for spin-orbit motion. In Proceedings to 
EPAC02, Paris and DESY-M-02-01 (2002) 

A.W. Chao: Evaluation of radiative spin polarization in [19 
an electron storage ring. Nucl. Instr. Meth. 180:29 (1981) 
Yu. Eidelman and V. Yakimenko: The application of Lie [20 
methods to the spin motion in nonlinear collider fields. 
Particle Accelerators, 45:17-35 (1994) 

S.R. Mane: Electron-spin polarization in high-energy [21 
storage rings. II. Evaluation of the equilibrium polariza- 
tion. Phys. Rev. A(36):120-130 (1987) [22 
K. Yokoya: Non-perturbative calculation of equilib- 



rium polarization of stored electron beams. KEK-92-6, 
Tsukuba (1992) 

V. Balandin and N. Golubeva: Nonlinear spin dynam- 
ics. Proceedings of the XV International Conference on 
High-energy Particle Accelerators, Hamburg, p. 998-1000 
(1992) 

K. Heinemann and G.H. Hoffstaetter: A tracking algo- 
rithm for the stable spin polarization field in storage rings 
using stroboscopic averaging. Phys. Rev. E 54:4240-4255 
(1996) 

G.H. Hoffstaetter, M. Vogt, and D.P. Barber: 
Higher-order effects in polarized proton dynamics. 
Phys. Rev. ST-AB, 2(11):1-15 (1999) 
D.P. Barber, G.H. Hoffstaetter, and M. Vogt: Polarized 
protons in HERA. In Proceedings to EPAC02, Paris and 
DESY-M-02-01 (2002) 

M. Froissart and R. Stora: Depolarisation d'un fais- 
ceau de protons polarises dans un synchrotron. Nucl. In- 
str. Meth. 7:297-305 (1960) 

V. Balandin, N. Golubeva, and D.P. Barber: Studies of 
the behavior of proton spin motion in HERA-p at high 
energies. DESY-M-96-04 (1996) 

Ya.S. Derbenev and A.M. Kondratenko: Acceleration 
of polarized particles. Sov. Phys. Doklady, 20:562, 1976. 
alsoinruss.: Dokl.Akad.Nauk Ser.Fiz. 223:830-833 (1975) 
SY. Lee: Spin dynamics and Snakes in synchrotrons. 
World Scientific (1997) 

D.P. Barber and A. De Roeck, editors: Proceedings of 
Polarized Protons at High Energies, DESY-Proceedings- 
1999-03 (1999) 

V. Anferov: 820 GeV/c polarized protons in HERA, 

SPINQHERA report update. In H|. 

A.W. Chao: Workshop summary - accelerator issues. In 

M. 



14 




100 300 500 700 

p [GeV/c) 




100 300 300 700 900 

p (GeV/c) 



FIG. 16: Improvement of the higher-order Pu m and v(J y ) by 
matching 8 snake angles and the orbital phases. The snake ar- 
rangement is (-|0000000)o (blue) and (^abc0-c-b-a)o (green). 
As a comparison Pu m from SODOM II is shown for the snake 
and phase advance matched four-snake scheme 0^^^)o (red 
background curve) . The resonances v — 2Q y and v = 1 — 2Q y 
are indicated (also red) 



[23] V. Ptitsin and Yu.M. Shatunov: The investigation of 
spin resonances in an accelerator with Siberian Snakes. In 
C.W. de Jager, T.J. Ketel, P.J. Mulders, J.E. J. Oberski, 
and M. Oskam-Tamoezer, editors, Proceedings SPIN96, 
p. 516-518. World Scientific (1996) 

[24] E.D. Courant: Hybrid helical snakes and rotators for 
RHIC. BNL-61920, Brookhaven (1995) 

[25] G.H. Hoffstaetter: Polarized protons in HERA. In 
DESY-96-05 (1996) 

[26] B.W. Montague: Polarized beams in high-energy storage 
rings. Physics Reports, 113(l):l-96 (1984) 

[27] G.H. Hoffstaetter: Matching of Siberian Snakes. AIP 
Conf.Proc.667:93-102 (2003) 



[28] G.H. Hoffstaetter: Snake matching. In UM-HE 99-05, 
University of Michigan (1999) 

[29] H. Rose: Hamiltonian magnetic optics. Nucl. In- 
str. Meth. A258:374-401 (1987) 

[30] R. Degenhardt: Korrektur von Aberrationen in der 
Teilchenoptik mit Hilfe von Symmetrien. Doctoral thesis, 
Darmstadt University of Technology (1992) 

[31] W. Wan: A Theory of Arbitrary Order Achromats. PhD 
thesis, Michigan State University (1995) 

[32] G.H. Hoffstaetter: Comments on aberration correc- 
tion in symmetric imaging energy filters. Nucl. In- 
str. Meth. A(427):275-281 (1998) 

[33] G.H. Hoffstaetter: Polarized protons in HERA. In Pro- 
ceedings of NUCLEON99, Nucl. Phys. A (1999) 

[34] D.P. Barber, M. Vogt, and G.H. Hoffstaetter: The 
amplitude-dependent spin tune and the invariant spin 
field in high-energy proton accelerators. In Proceedings 
EPAC98, Stockholm (1998) 

[35] G.H. Hoffstaetter: Successive approximations for 
charged particle motion. Ultramicroscopy 81 p. 111-121 
(2000) 

[36] Ch. Weifibacker and G. H. Hoffstaetter: Nonlinear spin 
transfer maps. In 120^1 . 

[37] K. Yokoya: An algorithm for calculating the spin tune 
in circular accelerators. DESY-99-006 (1999) 

[38] G.H. Hoffstaetter and M. Vogt: Sprint users guide and 
reference manual. DESY (2002) 

[39] D.P. Barber, K. Heinemann, G.H. Hoffstaetter, and 
M. Vogt: Calculations of the equilibrium spin distribu- 
tion for protons at high energy in HERA. In Proceedings 
of SPIN96, Amsterdam (1996) 

[40] D.P. Barber, G.H. Hoffstaetter, and M. Vogt: Using the 
amplitude dependent spin tune to study high order spin- 
orbit resonances in storage rings. In Proceedings SPIN00, 
Osaka (2000) 

[41] D.P. Barber, G.H. Hoffstaetter, and M. Vogt: Proton 
and electron polarsation in storage rings: some basic con- 
cepts. In Proceedings SPIN00, Osaka (2000) 

[42] D.P. Barber, G.H. Hoffstaetter, and M. Vogt: Numerical 
and theoretical spin stability studies for HERA-p. In |2oTj. 

[43] G.H. Hoffstaetter and M. Vogt: Odd and even order spin- 
orbit resonances with Siberian Snakes. In preparation for 
Phys. Rev. ST-AB (2004) 

[44] P. Wesolowski, K. Balewski, R. Brinkmann, Ya. Der- 
benev, and K. Flottmann: An injector study for electron 
cooling at PETRA using a bunched beam. In Proceedings 
of ECOOL99, Uppsala (1999) 

[45] M. Gentner, D. Husmann, C. Steier, R. Brinkmann, and 
Ya Derbenev: On the possibilities of electron cooling for 
HERA. Nucl. Instr. Meth. A424:277-295 (1999) 

[46] R. Brinkmann for the DESY electron cooling study 
group: Electron cooling studies for PETRA and HERA. 
In 

[47] D.P. Barber, G.H. Hoffstaetter, and M. Vogt: Spin mo- 
tion at high energy in HERA-p. In UM-HE 99-05, Uni- 
versity of Michigan (1999) 



